Specific heat and thermodynamic critical field for the molecular metallic hydrogen 
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In the framework of the Eliashberg formalism the free energy difference between the supercon- 
ducting and normal state for the molecular metallic hydrogen was calculated. The pressure values 
Pi = 347 GPa and P2 = 428 GPa were taken into consideration. It has been shown, that together 
£\) . with the increase of the pressure, grows the value of the specific heat jump at the critical temperature 

and the value of the thermodynamic critical field near zero Kelvin: [AC {To)] p2 I [AC (Tb)] pl — 2.33 
and [H c (0)] 2 / [Ha (0)] pl ~ 1.74. Next, it has been stated, that the ratio AC (To) /C N (To) also 
increases from 1.91 to 2.39; whereas T C C N (To) /-He (0) decreases from 0.152 to 0.140. The last 
results prove that the considered parameters significantly diverge from the prediction based on the 
BCS model. 
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The study of the metallic hydrogen's properties has been lasting for over seventy years. In 1935, Wigner and 
Huntington for the first time suggested, that under the influence of high pressure (p) the hydrogen should transform 
into the molecular metallic phase [Ik The later theoretical results set the metallization of hydrogen in the pressures 
range from 300 GPa to 400 GPa 0,0. It is worth mentioning, that the understanding of the high-pressure properties 
of hydrogen seems to be substantial due to the fact, that this element in the metallic state (both molecular and atomic) 
is appearing inside the planets of the Jovian type . 

The next step was made by Ashcroft who suggested, that the metallic hydrogen could be potential high-temperature 
superconductor [5]. Since that moment, the constant interest in the properties of the hydrogen's superconducting 
state has been dated. In particular, the numerical results predict that in the range of the "lower" pressures (up to 
500 GPa) the critical temperature (Tc) is of the order (80-300) K 0, 0, 0; 0- For the extremely hi gh p ressure 
(2000 GPa) the superconducting state in the atomic metallic hydrogen has been studied in the papers [jjj]. [lOj. It has 
been shown that the critical temperature decreases from 631 K to 413 K for \x* c 6 (0.1,0.5), where fi c denotes the 
■ critical value of the Coulomb pseudopotential. In the considered case the other thermodynamic parameters diverge 
I/-} \ from the BCS values [ll[ e.g.: the dimensionless ratio r\ = AC (Tc) /C N (Tc) is changing from 1.82 to 1.68 together 
lO ■ with the Coulomb pseudopotential's growth, whereas the minimum value of = TcC N (Tc) / H c (0) is equal to 0.162 
[Toj . The symbols defining the ratios r% and r-i have following meaning: AC (Tc) denotes the specific heat difference 
between the superconducting and normal state at the critical temperature, C N (Tc) represents the specific heat of 
the normal state, while He (0) is the value of the thermodynamic critical field at the temperature of zero Kelvin. 

In the literature the specific heat and the thermodynamic critical field were not determined for the molecular 
metallic hydrogen. Due to the large values of the electron-phonon constant ([A] = 0.93 and [A] = 1.2) it has to 
be presumed, that above quantities should be calculated in the framework of the Eliashberg formalism [T^. In the 
paper, we take into consideration the following values of the pressure: p\ = 347 GPa and p2 = 428 GPa. In this case 
the molecular metallic hydrogen crystallizes in the Cmca structure 0, [13{ . 



II. THE ELIASHBERG EQUATIONS 

The BCS theory is based on the Hamiltonian, which models the pairing interaction in the simplest effective way. 
We notice that the BCS Hamiltonian can be derived from the more realistic Frohlich's operator (Hp), which describes 
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the electron-phonon coupling in the open form llJl , [15| . The Eliashberg equations are derived directly from Hp with 
an use of the thermodynamic Green functions |l6j . As a result one can obtain 
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The solutions of the Eliashberg equations are two functions defined on the imaginary axis: the wave function renor- 
malization factor (Z n = Z (iu n )) and the order parameter function (cf> n = <f> (iw„)); w n = (ir//3) (2n — 1) is the n-th 
Matsubara frequency, where /3 = (fc^T) 1 (fee denotes the Boltzmann constant). In the framework of the Eliashberg 
formalism the order parameter is defined as: A„ = <f> n /Z n . The symbol A (z) represents the pairing kernel: 

A(z)^2 f W <m-^—ta 2 F(n). (3) 

The Eliashberg functions for the pressures pi and pi (a 2 F (f2)) were determined in the paper 0. The symbol Q max 
denotes the maximum phonon frequency, where [ri m ax] pi = 477 meV and [rimax] p2 = ^08 meV. 

The depairing correlations, appearing between electrons, are modeled with the help of the Coulomb pseudopotcntial 
fi^; the symbol 9 denotes the Heaviside unit function and w c is the cut-off frequency {uj c = 30 majc ). In the paper we 
have assumed low value of the Coulomb pseudopotential for both considered pressures ([i* c = 0.1). The assumption 
above can be justified by referring to the Bennemann-Garland formula [13]: Mc ~ 0-26p(0) / [1 + p(0)], where the 
symbol p (0) indicates the value of the electronic density of states at the Fermi energy. In particular, we have: 
pi (0) = 0.4512 states/Ry/spin for pi and p2 (0) = 0.4885 states/Ry/spin for p 2 @]- Thus, [Mclpi an( i [Mc] P 2 amoun ts 
- 0.081 and - 0.085 respectively. 

From the mathematical point of view the Eliashberg set is composed of the strongly non-linear algebraic equations 
with the integral kernel A (z) . In order to achieve stable solutions one needs to take into account adequately large 
number of the equations. In the paper we have assumed M = 800, what assured stability of the solutions beginning 
from the temperature of Tn = 11.6 K (1 meV). The Eliashberg equations were solved by using the iterative method 
presented in the papers [l8| and [l9j . 

III. THE NUMERICAL RESULTS 

The solutions of the Eliashberg equations for the selected temperatures have been presented in Figs. [T]and[2j It can 
be easily noticed, that the functions Z m and A m decrease together with the Matsubara frequencies' growth. However, 
Z m saturates considerably slower than A m . 

The applied pressure significantly influences on the values of wave function renormalization factor and the order 
parameter. From the physical point of view the above fact means, that together with the increasing of p increases the 
electron effective mass (m* ~ Z m= i) and the value of critical temperature (pc] = 108.2 K, [Tb] = 162.7 K). 

Analyzing the dependence of Z m and A m on temperature it has been stated, that the solutions of the Eliashberg 
equations very unlikely evolve with T. In Fig. [3] we have plotted the functions Z m= \ (T) and A m= i (T). The 
presented results show, that the wave function renormalization factor is weakly dependent on the temperature and 
takes its maximum for T = Tq. In contrast, the temperature dependence of the order parameter is strong and 

can be modeled by using the formula: A m= i (T) — A m= i (T ) — (^) > where: [A m= i (T )] pi = 18.15 meV, 

[A m=i (T )] p2 = 29.12 meV, [0\ pi = 3.58 and [/3] p2 = 3.61. 

The thermodynamic properties of the molecular metallic hydrogen can be exp licitly determined on the basis of the 
free energy difference between the superconducting and normal state (A.F) [20|: 
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FIG. 1: The wave function renormalization factor on the imaginary axis for selected values of the temperature. The figure (A) 
shows results for pi , the figure (B) for P2 . 




FIG. 2: The order parameter on the imaginary axis for selected values of the temperature. The figure (A) shows results for pi, 
the figure (B) for p2. 



where the functions and denote the wave function renormalization factors for the superconducting (S) and 
normal (N) state respectively. 

In the first step, on the basis of Eq. ((J]), we have calculated the specific heat difference between the superconducting 
and normal state (AC = C s - C N ) : 

AC = l rf 2 [AF/p(0)] 
W(0) P d(k B Tf ' 1 ' 

Next, the specific heat in the normal state has been calculated with an use of the formula: 

C N _ 7 
k B p(0) ~ 



(6) 



where 7 = |7r 2 (1 + A). In Fig. Uwe have plotted the temperature dependence of the specific heat for the super- 
conducting and normal state. Assuming previously given values of the electronic density of states it can be shown, 
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FIG. 3: (A) The dependence of the wave function renormalization factor for the first Matsubara frequency on the temperature. 
(B) The dependence of the order parameter for the first Matsubara frequency on the temperature. In both cases the results 
for p\ and p2 are presented. 

that together with the growth of p the specific heat's jump at the critical temperature very strongly increases. In 
particular, we have: [AC {T c )] p2 / [AC {T c )] pi ~ 2.33. 
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FIG. 4: The dependence of the specific heat in the superconducting and normal state on the temperature. The figure (A) 
shows results for pi, the figure (B) for p2- The vertical line indicates a position of the specific heat jump that occurs at Tc- 



Below, we have calculated the values of the thermodynamic critical field (cgs units): 



Vp¥) 



= y/-8n [AF/p(0)]. 



(7) 



In Fig. [5]we have presented the dependence of He I ' \J p (0) on the temperature. On the basis of obtained results we 
can see, that the value of the thermodynamic critical field near the temperature of zero Kelvin (He (0) ~ He (To)) 
also strongly increases with the pressure: [He (0)] P2 / [He (0)] pi — 1-74. 

On the basis of determined thermodynamic functions one can calculate two fundamental ratios: r\ and T2- Let 
us notice, that in the framework of BCS model these quantities have the universal values ([?~i] BC g = 1-43 and 
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FIG. 5: The thermodynamic critical field as a function of the temperature. The figure (A) shows results for pi, the figure (B) 
for P2 . 



[r2] BCS = 0.168) [TT(. For the molecular metallic hydrogen following results were obtained: 

[ri] pi = 1.91, N P2 - 2.39 (8) 

and 



[ra] pi = 0.152, [r 2 ] pa = 0.140. (9) 

It is easy to notice that the calculated ratios significantly diverge from the values predicted by the BCS theory. 
Additionally it should be underlined, that n is increasing together with the pressure's growth, whereas the ratio r 2 
is decreasing. 



IV. SUMMARY 



In the paper the free energy difference between the superconducting and normal state for the molecular metallic 
hydrogen was calculated. The pressure values p\ = 347 GPa and P2 — 428 GPa were taken into consideration. 
On the basis of achieved results it has been shown, that the specific heat's jump at the critical temperature and the 
thermodynamic critical field near the temperature of zero Kelvin strongly increase together with the pressure's growth 
([AC (T c )] p2 I [AC (T c )] pl ~ 2.33 and [H c (0)] p2 / [H c (0)] pl ~ 1.74). The obtained thermodynamic quantities enable 
the determination of the fundamental ratios: r\ and r 2 . It has been proven, that the ratios t\ and r 2 very considerably 
differ from the values predicted by the BCS model. In particular, T\ is increasing from 1.91 to 2.39 together with the 
pressure's growth; whereas r 2 is decreasing from 0.152 to 0.140. 
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